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Abstract 

Thermodynamical fluctuations of temperature in mirrors of gravitational 
wave antennae are transformed through thermal expansion coefficient into 
additional noise. This source of noise, which may also be interpreted as fluc- 
tuations due to thermoelastic damping, may not be neglected and leads to 
the necessity to reexamine the choice of materials for the mirrors. Additional 
source of noise are fluctuations of the mirrors' surfaces caused by optical power 
absorbed in dielectrical reflective layers. 

I. INTRODUCTION 

The full scale terrestrial gravitational wave antennae are now in the process of assembling 
and tuning. It is expected that in the year 2001 in one of these antennae (project LIGO- 
I) the sensitivity in units of amplitude of the perturbation of metric will be at the level 
/i ~ 1 X 10^^^ IIIJH- This value means that the optical readout system will be able to register 
small displacement between two test masses (two heavy optical mirrors placed in vacuum 
pipe at the distance L = 4 x 10^ cm) with the resolution ALgrav — ^ hL 2x 10^^^ cm. This 
displacement may be produced by a burst of gravitational wave generated in astrophysical 
catastrophe. The characteristic frequency in this wave packet and, correspondingly, in the 
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displacement cjgrav has to be of the order of 27r x 10^ s~^. The improvement of the readout 
system and better isolation of the mirrors from different sources of noise will permit to 
achieve higher sensitivity: in the year 2005 the planned level of sensitivity has to be better 
than ^ 2 X 10~^^ and correspondingly it will be necessary to register ALgrav < 4 x 10^^^ cm. 
It is worth noting that in the prototype LIGO (where L = 4 x 10'^ cm) the already achieved 
resolution was ALgrav = 3 x 10~^^ cm f^. 

The important feature of this antenna is that the optical readout system does not register 
displacements between the centers of the mirrors' masses but between the internal surfaces of 
the mirrors. These surfaces are coated by several thin dielectrical layers, which provide high 
optical reflectivity of the mirrors and high sensitivity to the displacements. The necessity to 
measure the displacements between the internal surfaces inevitably determines the require- 
ment for high quality factors Qint (small loss angle) of internal mechanical modes of the 
mirrors, which provide relatively small values of random thermal fluctuations of the mirrors 
related to the mirrors' center of mass because the value of Wgrav is substantially smaller than 
the eigen-frequencies of internal modes. We shall call as Brownian the fluctuations which are 
usually calculated from phenomenological model of structural losses in the material in order 
to distinguish them from another type of thermodynamical fluctuations which we describe 
below. This requirement for high Qmt determines the choice of the mirrors' material. On 
the first stage of LIGO the mirrors are made of very pure fused silica {Si02), because the 
values of Qintsi02 exceed 2 x 10^ Oil- 

At present the possibility to make the mirrors for the second stage of LIGO from sapphire 
monocrystal (^4/203) is widely debated, because more than 20 years ago the achieved value 
of QintAbOs for this material was 3 x 10^ [|]. 

We show below that the value of Qint is not the only parameter which determines the 
choice of material from which the mirror has to be manufactured. Till now the experimental- 
ists' interest has neen concentrated on the reduction of the Brownian motion of the mirrors' 
surfaces (as it was described above), to the reduction of thermal fluctuations of the centers 
of mirrors' masses due to the friction in the mirrors' suspension and to excess noise in 
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the suspensions To our best knowledge nobody has so far paid attention to nonhnear 
effects which may be "responsible" for additional fluctuations of the mirrors surface, which 
may exceed traditional Brownian noise for low frequencies. The simplest of these effects 
is thermal expansion, originated from anharmonisity of interaction potential of atoms in 
solid state. The following numerical example illustrates the role of nonzero value of thermal 
expansion coefficient a of the mirrors material. Suppose that the surface layer of the mirror 
having thickness It — 10"^ cm (it is the characteristic length of the decay of Fourier thermal 
wave at frequency 27r x 10^s~^ in fused silica) is heated due to some process on 6T. Then 
the distance between the surface of the mirror and its center of mass will also change on the 
value of 61: 

5l = al5Tc^lx 10-^^ cm X f -\ x { \ x { ^^^^^ . (1) 

V5xlO-7K-V \lxl0-2cmy V2xlO-9Ky 

The estimate (|1]) shows, that if the mirrors are made of fused silica and if the temperature 
of the first layer (10~^ cm thick) of the mirror changes by 2 nK during 5 x 10^'^ s then the 
antenna will register a "burst" of gravitational radiation with amplitude of the order of 
h ~ 10^^^ with cUgrav = 27r X 10^ s. 

The aim of this article is to present the results of the analysis of the two effects which 
are "responsible" for random fluctuations of temperature in the bulk of the mirrors and 
which may mimic at certain level the gravitational wave if a 7^ 0: photo-thermal shot 
noise (due to random absorption of optical photons in the surface layer of the mirror) and 
thermodynamical fluctuations of temperature. 

As in both effects we are interested in the fluctuations of the coordinate averaged over 
the spot on mirror surface with radius rg ^ 1.5 cm of laser beam, which is much smaller 
than the radius of mirror ~ 15 cm, we replace the mirror by half-space: < x < 00, 00 < 
y < 00, cx3<z<c)0. We consider also for simplicity this half-space to be isotropic (the 
anisotropy of different constants for sapphire is not very large). The role of these two effects 
for the condition to reach the standard quantum limit of frequency stability was analyzed 



more than 20 years ago [|T0|. 



As at room temperature the power radiated from the surface (Stephan-Boltzmann law) 
is much lower than the heat exchange due to thermal conductivity (Fourier law), below we 
use the simplified boundary condition: 

du{x,y,z,t) ^ ^ 

x=0 

where u is the deviation of temperature from the mean value T. 

For the calculations we use approximation in which we neglect the effects of limited speed 
of sound but take into account thermal relaxation. 



II. PHOTO-THERMAL SHOT NOISE IN THE MIRRORS OF THE ANTENNAE 

The multilayer coating of the mirrors absorbs a fraction of power W of the optical beam, 
circulating between the mirrors of Fabry- Perot resonator. This fraction for the best existing 
today coatings has not been measured with good precision yet, but for estimates the usually 
accepted value of this fraction is approximately equal to 5 x 10^'' ^ 10^^ [|11|]. The value 



of W sufficient to reach h ~ 10^^^ (which is close to Standard Quantum Limit (SQL) of 
sensitivity) is approximately 14^ ~ 1 x lO^Watt = lO^^erg/sec or A^^o — 5 x 10^^ photon/s 
|T^. It is reasonable to estimate that each mirror absorbs Wabs — 1 Watt = 10^ erg/s or 



A^ ~ 5 X 10^^ photon/s. Each absorbed optical photon having energy huo ~ 2 x 10^^^ erg 
gives birth to a bunch of approximately 50 thermal phonons. Each of these phonons has a 
relatively short free path /* (for Si02 and AI2O3 the values are lsi02 — 8 ^ 10^^ cm and 
l^i^Os — 5 X 10^^ cm ||13[). After a very short time interval which is a few times longer than 
l*/vs — 10^^^ s (where Vg is the speed of sound) these phonons will produce a local jump of 
temperature inside the coating (5 -^ 10) x 10^'^ cm thick. As we want to know the variance 
of temperature during the averaging time vr/cUgrav ~ 5 x 10^'^ s which is several orders longer 
than l*/vs and as during this time a large number (A^vr/tUgrav ~ 2.5 x 10^^) of photons are 
creating bunches of thermal phonons — it is reasonable to use the model of shot noise. As 
the length 1^ is much smaller than the radius of laser beam tq ~ 1.5 cm for semi-qualitative 
analysis one can use one dimensional model, described by thermal conductivity equation: 



du{x, t) 



dx 



= 0, < X < oo, (4) 

z=0 

wit) wit')) = hu^W^Sit - t'), wit) = iy,bs - Wo, Wo = {WaUt)), 



where = A*/ (pC), A* is thermal conductivity, p is density and C is specific heat capacity. 

Assuming that displacement Xid (which is of interest) is proportional to temperature u 
averaged along the axis x 

roo 

= a / dx uix,t). (5) 

JO 

it is not difficult to obtain the spectral density of Xid Q: 

c / N o 2 ^^qWo 1 , 
5'tsidH =2a — — — 2-^— . (6) 



For the exact calculations for half-space (see Appendix A) we introduce the displacement 
X averaged over the beam spot on the surface 

^ = ^1 r dydz,v^ix = Q,y,z)e-^y"^-^"^''l (7) 
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Trrg J J—oo 



where v^ix = 0,y,z) is x-component of vector of deformation on the surface, and find the 
spectral density of X: 



ipC-KTo)^ ijj^ 



where a is the Poisson coefficient. 



^ We use "one-sided" spectral density, defined only for positive frequencies, which may be calcu- 
lated from correlation function {Xit)Xit -|- r)) using formula 

/oo 
dr(X(t)X(t + r)) cos(cjr). 
-co 
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Here and below we use the fact that for materials at room temperature the thermal 
relaxation time of the spot tt ^ rg/a^ is very large: 

u » 3. (9) 

^0 

It is easy to see that the exact solution differs from the approximate one (|^) only by 
Poisson coefficient. 

III. THERMODYNAMICAL FLUCTUATIONS OF TEMPERATURE, 
THERMOELASTIC DAMPING AND SURFACE FLUCTUATIONS 



It is known from thermodynamics that the total variance of fluctuations of temper- 
ature in volume V is described by the following simple formula: 

^^T^)-^. (10) 

where n is the Boltzmann constant. In classical solid state thermodynamics these fluctua- 
tions are not correlated with fluctuations of volume (which are responsible for the Brownian 
noise) |]14[ but this is not the case for nonzero thermal expansion coefficient. To find the 
influence of the fraction of this type of fluctuations in the vicinity of Ug^av on the vibration 
of the surface along the x-axis in direct way we use the Langevin approach and introduce 
fluctuational thermal sources F{r,t) added to the right part of the equation of thermal 
conductivity: 

^-a'Au = F{r,t). (11) 

These sources should be normalized in a way to satisfy formula ([T0|). By substituting the 
solution of this equation for the temperature u into the equation of elasticity and averaging 
displacement (0) over the spot on the surface, we find spectral density of its fluctuations in 
frequency range of interest (see details in Appendix B): 



We can also find the spectral density of displacement using Fluctuation-Dissipation the- 
orem (FDT) if we assume that the only dissipation mechanism in the mirror is thermoelastic 
damping. Elastic deformations in a solid body through thermal expansion lead to inhomo- 



geneous distribution of temperature and hence to fiuxes of heat and losses of energy [19 



The same approach was used IT^JTH] for the calculation of vibration of the surface due to 



the Brownian fiuctuations. To calculate the fiuctuations in this approach we should apply 
a periodic pressure p distributed over the beam spot on the surface: 



z,t) = ^ e-^J^'+^'^/^o e^-* (13) 

and we should find susceptibility x = X /Fq. After that in accordance with FDT the spectral 
density of displacement X must be proportional to imaginary part of susceptibility x- 

5.H = ^|Im(xH)|. (14) 

UJ 



Calculations (see Appendix C) show that FDT gives the same result as formula (|T2|). 
Therefore, we can conclude that thermodynamical fiuctuations of temperature are the phys- 
ical source of fiuctuations deduced from FDT based on thermoelastic damping. 



Calculations of Im(x) are rather bulky but Zener suggested a very simple approach 
allowing to estimate thermoelastic loss angles and hence the imaginary parts of elastic co- 
efficients. It allows to find Im(x) only from formula for x in zero order approximation (see 
Appendix C): 

= <i^, (15) 

where E is the Young's modulus. 

The imaginary parts of elastic coefficients due to thermoelastic damping may be calcu- 



lated in Zener 's approach according to the formulas [|l^]: 

Im(cr) as — ctt TE 1 + a uott 



a ax 1 + cj^Tj. Cp a 1 + u'^Tj,^ 

lm{E) Es — Et ujtt o? TE ujtt 



E Et 1 + u;V| Cp 1 + cuVf 
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(16) 



Where Et, Es are consequently iso-thermic and adiabatic Young's moduluses {Et ^ Es — 
E) and axiC^s are iso-thermic and adiabatic Poisson coefficients {ut ^ cts — cr), tt is the 
time of thermal relaxation of our spot. By substituting these values in (|I5|) and then in 
(0) and taking into account that tt = 2rQ/a^ ^ l/cj (this relaxation time may be found 
relatively easy from the solution of the problem of thermal relaxation of the heated spot) 
we finally obtain the same formula (p!2D . 

It is really worth noting that all the three methods give precisely the same result. 



IV. NUMERICAL ESTIMATES 

Now we want to compare the fluctuations of the mirrors' surface caused by the described 
above mechanisms with other known sources of noise in LIGO antennae. 

The main determining source of noise associated with the mirrors' material is usually 
described by the losses in the model of structural damping (we denote it as Brownian motion 
of the surface). In this model the angle of losses (p does not depend on frequency and 
analogously to P, p!7| , p^ by substituting E = E{1 + icp) in (p!5| ) and then in (|I^ we obtain: 



S3oM^^^<^^. (17) 



o 



For simplicity here we neglected the unknown imaginary part of Poisson coefficient. 

The sensitivity of gravitational wave antenna to the perturbation of metric may be 
recalculated from noise spectral density of displacement X using the following formula: 

M-) = (18) 

where we used the fact that antenna has two arms (with length L) and in each arm only 
one mirror adds to noise (on end mirrors the fluctuations are averaged over the larger beam 
spot). 

The LIGO-II antenna will reach the level of SQL, so we also compare the noise limited 
sensitivity to this limit in spectral form ||12[|: 



8 



On figures 1 and 2 we plot the sensitivity limitations for the perturbation of metric 
caused by the described thermal mechanisms and SQL (formulas |],|12|,|17|,|18|,|19]) for fused 
silica (Fig.l) and sapphire (Fig.2) mirrors. We used the following numerical values for the 
parameters of the two materials — fused silica and sapphire. (We should note that the actual 
figures, which various sapphire manufactures provide, significantly vary - tens of percents, 
especially for the thermal expansion coefficient): 

= 27r X 100 s-\ ro = 1.56 cm, T = 300 K, 

Wo = 2 X 10^^ s-\ L = 4:X 10^ cm, Wq = lOVg/s; (20) 

Fused silica: a = 5.5 x 10"^ K-\ X* = 1.4 x 10^ (21) 

cm s K 

p = 2.2-^, C = 6.7x10^-^, m = l.lxl0^g, 
cm'^ g K 

E = 7.2x10"—^, a = 0.17, = 5x10"*^; 
cm'^ 

Sapphire: a = 5.0 x 10"^ K-\ X* = 4.0 x 10^ , (22) 

cm s K 

p = 4.0 C = 7.9 X 10^ m = 3 X 10^ g, 

cm'^ g K 

E = Ax 10^2 a = 0.29, = 3 x 10"^ (23) 

cm'^ 

From these two figures it is evident that the effects associated with thermal expansion in 
no case may be neglected. Especially unfortunate result is obtained for sapphire - the noise 
from thermodynamical fluctuations of temperature dominates for frequencies of interest and 
is several times larger than the SQL. While for fused silica the situation is the opposite. The 
SQL is one order larger than TD fluctuations. It is important to note, that unlike structural 
damping for which angles of losses at low frequencies were not yet directly measured, our 
mechanism has fundamental nature and is calculated explicitely without model assumptions. 
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CONCLUSION 



The analysis and numerical estimates presented above for two effects allow to give some 
recommendations for the strategy of upgrading laser interferometric gravitational wave an- 
tennae. 

1. It is important to emphasize that contrastingly to the gravitational wave, both effects 
do not displace the centers of masses of the mirrors. Thus, in principle it is possible to 
subtract these effects as well as usual Brownian fluctuations (with certain level of accuracy) 
in a way similar to the procedure suggested for the subtraction of thermal noise in suspen- 
sion fibers However, till now to our knowledge nobody presented the scheme of such 
subtraction. 

2. Potentially there exists a possibility to decrease substantially these effects by choos- 
ing materials for the mirrors with substantially smaller values of a. There also exists the 
possibility to use special cuts of anisotropic monocrystals for which a turns to zero. 

3. Thermodynamic fluctuations decreese with growing radius of the beam spot as r^J"^ 
even stronger then for the Brownian fluctuations. In this way, some optimization of the 
geometry of the resonators may help. 

4. The "brute force" method for the improvement of the sensitivity by increasing the 
circulating optical power W does not look promising, taking into account photo-thermal 
shot noise. It seems that few mega- Watts is the upper limit for W . 



ADDENDUM 



The evident method to suppress the thermodynamical fluctuations by using larger 
beamspots is unfortunately limited by diffraction losses and technology. If we choose the 
LIGO-I geometry with beam sizes on internal and end mirrors: r^i = 3.6/^/2 cm, roe = 
4.6/-\/2 cm then the threshold of sensitivity of thermodynamical fluctuations will be only 
1.7 times lower than on our graphics still 2 times higher than the SQL. 
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APPENDIX A. 



In this Appendix we calculate the spectral density (|]) of the fluctuations of the averaged 
displacement X caused by photo-thermal shot noise. 
We start from thermal conductivity equation: 



— a Au ~" 

dt 



pC tcvq 



2(5 (x) exp 



(24) 



{Wit)wit')) = hiUoWodit - t'), wit) = lyabs(t) - Wo, Wo = (W^abs), 



with boundary condition (0) and equation of elasticity for the field of deformation v{x, y, z) 
(see [p!5tp^ ) with boundary conditions for tensions on free surface: 
1 



■ grad div v — — r-rot rot v = a grad m, 



cr 



E 



1 - 2a 
E 

2{l + <y) 



2(1 + a) 
(y ( dv^ _^ dvy ^ dv^ 
1 + a \dx dy dz 



au + 



2cr dv^ 



1 + a dx 



0, 



dy dx 



0, 



x=0 

E ( dvx dv 



x=0 



dz ~^ dx 



(25) 
(26) 
0, (27) 



x=0 



2(1 + a 

The problem (|2^ , |^) in half-space can be replaced by the problem ([2^ ) in full space, in 
which the function of source is evenly continued for negative x. In this case the condition (^ 
is satisfied automatically. Assuming that the absorbed power changes as w{t) = t/;(a;)e*'^* 
one can find the spectral component: 



°° dk^dkydkz ik^x+ikyy+ik^z „~r^(kl+kl)/4, 2w{uj) 



, (27r)3 pC{a^k^ + iuj)' 

After substitution of u{r,u!) into (^, ^B|) one can find the solution for -u as a sum ||T^ 
V = grad (p + v^"'\ where (p satisfies Poisson's equation (it can be deduced from (|25D ): 
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Ao9 = au 

1 — cr 



(28) 



without boundary conditions and function w^") satisfies equation (^) with zero right part 
and boundary conditions: 

E 



l + cr 



Qy2 QyZ 



(29) 



x=0 



a. 



E 9V 



xy 



(1 + a) dxdy 



E 9V 



x=0 



x=0 



(1 + (t) dxdz 

It is convenient to choose the function in symmetrical form (by evenly continuing the 
u{r,uj) for negative values of x): 

a{l + a) f f f°o dx'dy'dz' u{r,uj) 



(p{x,y,z,uj) 
ip{x = 0,y,z,uj) 



47r(l -a) J J ^(^x'y^{y-y'y + {z-z'y' 



(30) 



ail + a)2wiuj) f f dk^dkydk, ^_,.ii^kl+kl)/A+ikyy+ik,z 



[I -a) pC J J J-oo{2T:fk^{a?k^ + iuj) 
In this case dx<p{x = 0,y,z) = and therefore (i) the function gives no contribution into 



X and (ii) tangential tensions a^y and a^z in boundary conditions ( pQ]) turn to zero. 

The solution of the problem for v^"-^ is known (see for example [|15|)- Using condition 
and following table integrals: 



/ / dydz 

J J —oo 



v/p"+^ k±' 



k'y + kl 



yo ^ikxX+ikyy+ikzZ 

dxdydz^^ = jj, k= Jk^ + k^ + k^ 

oo y/x^ + y^ + z^ k'^ * 



(31) 
(32) 



we substitute v^""^ into (ffh and obtain: 



X(uj) 



-4a(l + a) 



W[LU) 



pC 

a(l + cr) w{u!) 



dJvxdkydJv^ 



g -r2(fc2+fc2)/2 



(27r)3 k^ia'^k^ + iu) 

" rffc,. /~ dA:^ e-^o'^i/^ = _^il±^ !f!M 1. (33) 

ivr^ci; pC i-oo ^ Jo kl + k\ in pCrl uo 

Assuming that spectral density of absorbed power is equal to Sw{uj) = 2huJoWo one can find 
the formula (BD- 
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APPENDIX B 



In this Appendix we deduce the formula (0) using Langevin abroach for thermal con- 
ductivity equation ([TT| ) and recalculate the thermo-dynamical fluctuations of temperature 
into deformations. 

Using time and space Fourier transform one can write down the formal solution of (fTJ): 

-(r,t) = ir§^^§^e~'*-', (34) 

We interest in temperature u averaged over some volume V. Assuming that 
( F{k, uj)F*{k', io') ) = (27r)4 F^\k\^6{k - k') 5{uj - u') one can calculate ( ) 
1 f ^^1 f f f f r dkdk'dujduj' {F{k,u)F*{k',u')) 



V Jv V Jv J J J J-oo (27r)8 {a^{ky + tuj){a^{k'y + tuj') 

Here we equate (pS]) to (|TOD and find Fq. 

In order to find TD temperature fluctuations in half-space with boundary condition (0) 
one can solve equivalent problem for full space with source function even continued for 
negative x, so that its correlator is equal to: 

{F{k,uj)F*{k',u') ) = {27cy F^\k\^ [6{k, - k',) + 5{k., + k',)] 6{ky - k'^) 5{K - K) 6{uj - u') 

After that we must substitute the expression (Q) for u into elasticity equation (|25| , 
pUI) . The solution for v can be found according to the same scheme as in Appendix A: 
V = grad if + v^"'\ where satisfies Poisson equation p8|). Again can be written in 
symmetrical form ( pO]) and the function v^^^ must satisfy equation p5| ) with zero right part 
and following boundary conditions: 

a^^ix = 0, y, z) = 1 1 1 r ^(^'^) K + 



[l-a) J J J J-oo (27r)4 a2A;2 + p 
c^xy(2; = 0, y, z) = a^z{x = 0, y, z) = 0. 



The solution for v^""^ is known W5\ 



= 0, 2/, = — / / dy'dz' 



nE 



—oo 



^[y-y'Y + ^Z-z'Y 
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Substituting it into (|^ we find X{t) and spectral density SToi^)- 

X{t) = 2a(l + a) (Ur dk.dykj^ JiM. % 
^ ' ^ ' n ] Loo (27r)4 a^k^ + iuj 

C 0^2.. , ^2 r f r dk^dkydK 2{kl) ,^,2/2 

Std[uj) - Sa [l + a) JJJ_^ ^2^^3 ^4^4 + ^2 

From the last formula, using integration in spherical coordinate system and neglecting the 
term a^k^ in denominator according to condition (P), we obtain formula (12) asymptotically 
correct for the frequencies ujr1/a? ^ 1. 



The table integrals ^2]) was used above in calculations of axx{x = 0,y,z), t>^")(x 



APPENDIX C 

To calculate the fluctuations of the surface using FDT approach with thermoelastic losses 
we should solve the system: 

1-(T,,^1-2(T 

grad div v — — -rot rot v = agrad u, (36) 

1 + 0" 2(1 + a) 

du 2 a aET d div v 

" " = Cp(l - 2a) dt ' 



^ XX xy xz ^ n 

tttq ox 



0. (38) 



x=0 



In zero order approximation we simply solve elasticity problem neglecting temperature u in 
the right part of the first equation in (|36D and in boarder conditions. This solution is well 



known ||15|| : 



Foil + (t) 



27r2£;rg 



00 



—00 

2' 



Girl = [-^^ + ^) +[-3- ^:^)) (y + ^ ^^)- (39) 

By substituting this solution in (^ and using (^Tj) we obtain susceptibility x^^'^ (0)- This 
result totally coincides with the calculations for the half-space in [^] (if we note that the 
stored static strain energy U = X(o)Fo/2 = F2x(°V2 and Wo = V2ro) and has minor 
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numerical difference with the approximation in |]T^. To find first order approximation for 
the distribution of temperature we substitute 

in (^) to find that 

(1) _ 4aT(l + a)Foiuj r r r-- dk, dky dk, ^-k^,rliA+ikr /..n 

Cp 7 7 7-00 ( 27r)3 A;2(a2A;2 + ,^)^ • ^^^^ 

We then substitute this solution in (^) to find second order approximation in the same 
way as in the two Appendices above by finding appropriate evenly continued ip function 
satisfying ( pSj |^). We find using (|32D 



^ A - 4(l+0-)a^TEFoZ^ /•/•/- dk, dky dk, kj k^,r^J,+^k,y+^k.z 

^ (l-a)pC 7 7 7-00 (27r)3 fc4(a2A;2 + ,,^)^ 



(42) 



By substituting this expression in 



(o,i/,2;) = — ^ dydz I =, (43) 



oo 



and using (^ we may finally obtain 



.... . _3,1 . ... ///_: ""'^'^J.^f ^'^ (44, 

and hence Im(x(u;)) which leads through (|T^ to precisely the same expression as (|T2p. 
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Sapphire (AI2O3) 
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Figure captions: 

Fig.l The sensitivity limitations for the perturbation of metric caused by the SQL and 
by thermal fluctuations in fused silica mirrors. 

Fig. 2 The sensitivity limitations for the perturbation of metric caused by the SQL and 
by thermal fluctuations in sapphire mirrors. 
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